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At high-energies the gluon-gluon fusion is the dominant mechanism of cc production. This 
process was calculated in the NLO collinear as well as in the k t -factorization approaches in the 
past. In this presentation we concentrate on production of cc pairs including several sublead- 
ing mechanisms. This includes: gg — > QQ, 7<? — > QQ, 91 — >• QQ, 77 QQ- In this context 
we use MRST-QED parton distributions which include photon as a parton in the proton as 
well as elastic photon distributions calculated in the equivalent photon approximation. We 
present distributions in the c quark (c antiquark) rapidity and transverse momenta and com- 
pare them to the dominant gluon-gluon fusion contribution. valent photon approximation. We 
discuss also inclusive single and central diffractive processes using diffractive parton distribu- 
tion found from the analysis of HERA diffractive data. As in the previous case we present 
distribution in c (c) rapidity and transverse momentum. Next we present results for exclusive 
central diffractive mechanism discussed recently in the literature. We show corresponding 
differential distributions and compare them with corresponding distributions for single and 
central diffractive components. Finally we discuss production of two pairs of cc within a 
simple formalism of double-parton scattering (DPS). Surprisingly very large cross sections, 
comparable to single-parton scattering (SPS) contribution, are predicted for LHC energies. 

1 Introduction 

In this presentation we discuss contributions of some subleading mechanisms neglected in the 
analysis of cc production. We include contributions of photon- gluon (gluon-photon) as well as 
purely electromagnetic photon-photon fusion. Here we present only some selective results. The 
formalism and more details has been shown and discussed elsewhere^. 

We discuss also diffractive processes (single and central) in the framework of Ingelman- 
Schlein model corr ected fo r absorption. Such a model was used in the estimation of several 
diffractive processes^EEIS 

The absorption corrections are necessary to understand a huge Regge-factorization breaking 
observed in single and central production at Tevatron. 



2 Production of heavy quarks 

The cross section for the cc production, assuming gluon-gluon fusion, was calculated both in 
collinear and k t fa ctori zation approaches. Our group has done detailed calculations in the second 
approach (see e.g. ^E). 

In the leading-order (LO) approximation within the /ct-factorization approach the quadruply 
differential cross section in the rapidity of Q (yi), in the rapidity of Q (2/2) and in the transverse 



momentum of Q (pi,t) and Q (j)2,t) can De written as' 

da / d 2 Kl,i d 2 K2, 



dyidy2d 2 pi t td 2 p2 t t ^— ' / w w 167r 2 (a:iiC2s) 2 



where Fi{xi,k{ t ) and J r j(x2, «| t) are the so-called unintegrated gluon (parton) distributions. 
The unintegrated parton distributions must be evaluated at: 

xi = -^exp(yi)H ^exp^J, 

^2 = — ^exp(-2/i)H ^exp(-y 2 ), 



where rrij t = \lv\t + m o 



3 Photon induced production of heavy quarks 

The dominant contributions of heavy quark-antiquark production are initiated by gluon-gluon 
fusion or quark-antiquark annihilation. In general, even photon can be a constituent of the 
proton. This idea was considered by Martin, Roberts, Stirling and Thorne in RefP. 

If the photon is a constituent of the nucleon then other mechanisms of cc production pre- 
sented in FigQ] are possible. 




Figure 1: Diagrams representing photon-induced mechanisms of heavy quark production. 



4 Results 

4-1 Gluon-gluon fusion 

Before we go to the new mechanisms we will present results for the dominant gluon-gluon fusion. 
In Fig|2] we show distributions in transverse momentum of c (or c) for the gluon-gluon fusion 
mechanism for different popular choices of scales (p? = 4m^, M^ 5 ,pf + m^). We show our results 
for ^/s = 500 GeV (left panel) and ^fs = 14 TeV (right panel). In this calculation we have used 
GRVS PDFs. The figure shows typical uncertainties due to the choice of the scale. We wish 
to stress here that at the higher energies the results of the calculations depend on the gluon 
distributions at small values of x. 



4-2 "fg and g^ subprocesses 



In Fig|3] we show transverse momentum distributions for the dominant gluon-gluon as well as 
for the subleading photon-gluon (gluon-photon) and photon-photon components for different 




Pt (GeV) p, (GeV) 



Figure 2: Distribution in quark/antiquark transverse momentum at y/s = 500 GeV (left panel) and for y/s = 14 
TeV (right panel) for different choices of scales and for GRV gluon distribution. 



gluon distribution functions EMEU f or the RHIC energy y/s = 500 GeV and for the nominal 
LHC energy y/s = 14 TeV, respectively. At the LHC energy the results for different GPDs differ 
considerably which is a consequence of the poorly known small-x region. The differences at the 
energy y/s = 14 TeV are particularly large which can be explained by the fact that a product of 
gluon distributions (both at small x) enters the cross section formula. New measurement of cc 
at the nominal LHC energy will be therefore a severe test of gluon distributions at small x and 
not too high factorization scales not tested so far. Similar uncertainties for the jg and 37 are 
smaller as here only one gluon distribution appears in the corresponding cross section formula. 
The uncertainties for the photon distributions are not yet quantified. 




Figure 3: Transverse momentum distribution for the standard gluon- gluon mixed gluon-photon and photon-gluon 
as well as for photon- photon contributions for RHIC (left) and LHC (right). 



It is very difficult to quantify uncertainties related to photon PDFs as only one set of PDFs 
includes photon as a parton of the proton. Here the isospin symmetry violation (not well known 
at present) would be an useful limitation. Our collection of the results for the photon induced 
mechanisms show that they are rather small and their identification would be rather difficult as 
the different distributions are very similar to those for the gluon-gluon fusion. Our intension is to 
document all the subleading terms. Our etimation shows that the sum of all the photon induced 
terms is less than 0.5 % and is by almost 2 orders of magnitude smaller than the uncertainties 
of the dominant leading-order gluon-gluon term. 



5 Single and central diffraction 



5. 1 Formalism 



The mechanisms of the diffractive production of heavy quarks (cc) are shown in FigsUl The 
formalism how to calculate respective cross section has been presented elsewhere^. 





Figure 4: The mechanism of single-diffractive production of cc. 




Figure 5: The mechanism of central-diffractive production of cc. 



5.2 Results 

In Figj6] we show transverse momentum distributions of charm quarks (or antiquarks). The 
distribution for single diffractive component is smaller than that for the inclusivegluon-gluon 
fusion by almost 2 orders of magnitude. Our results include gap survival factor M The cross 
section for central diffractive component is smaller by additional order of magnitude. 
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Figure 6: Transverse momentum distribution of c quarks (antiquarks) for RHIC energy y/s = 500 GeV (left panel) 
and for LHC energy y/s = 14 TeV (right panel) for the GRV94 gluon distributions. The result for single diffractive 
(Od or dO), central diffractive (dd) mechanisms are compared with that for the standard gluon-gluon fusion (00). 



In FigJT] we show distributions in quark (antiquark) rapidity. We show separately con- 
tributions of two different single-diffractive components (which give the same distributions in 



transverse momentum) and the contribution of central-diffractive component in Figj6j When 
added together the single-diffractive components produce a distribution in rapidity similar in 
shape to that for the standard inclusive case. 




Figure 7: Rapidity distribution of c quarks (antiquarks) for RHIC energy y/s = 500 GeV (left panel) and for 
LHC energy yfs = 14 TeV (right panel) for the GRV94 gluon distributions. The result for single diffractive (Od 
or dO), central diffractive (dd) mechanisms are compared with that for the standard gluon-gluon fusion (00). 



The cross section for single and central diffraction is rather small compared to the dominant 
gluon-gluon fusion component. However, a very specific final state should allow its identification 
by imposing special conditions on the one-side (single-diffractive process) and both-side (central 
diffractive process) rapidity gaps. We hope that such an analysis is possible at LHC. Special 
care should be devoted to the observation of the exclusive cc production. Without a special 
analysis of the final state multiplicity the exclusive cc production may look like an inclusive 
central diffraction. 

6 Production of two cc pairs in double-parton scattering 

The general formula for the cross section in terms of double-parton distributions can be written 

El 

da DPS = —^Fgg(x 1 ,x 2 ,f4,f4) F gg( x i x 2,f j2 i,f4) 

2(J eff 

da gg ^ c5 (x 1 ,x' 1 ,fil)dagg^ C c(x 2 ,x' 2 ,f4) dx 1 dx 2 dx' 1 dx' 2 . (1) 

In Fig. [H] we compare cross sections for the single and double-parton scattering as a function 
of proton-proton center-of-mass energy. At low energies the conventional single-parton scattering 
dominates. At low energy the cc or cccc cross sections are much smaller than the total cross 
section. At higher energies the contributions dangerously approach the expected total cross 
section^. This shows that inclusion of unitarity effect and/or saturation of parton distributions 
may be necessary. The effect of saturation in cc production has been included but not checked 
versus experimental data. Presence of double-parton scattering changes the situation. At LHC 
energies the cross section for both terms become comparable^. This is a completely new situation 
when the double-parton scattering gives a huge contribution to inclusive charm production. 

In Fig. [9j we present single c (c) distributions. Within approximations made in this paper 
the distributions are identical in shape to single-parton scattering distributions. This means that 
double-scattering contribution produces naturally an extra center-of-mass energy dependent K 

"New experiments at LHC will provide new input for parametrizations of the total cross section. 
6 If inclusive cross section for core was shown the cross section should be multiplied by a factor of two - two 
c or two c in each event. 




Figure 8: Total LO cross section for single-parton and double-parton scattering as a function of center-of-mass 
energy (left panel) and uncertainties due to the choice of (factorization, renormalization) scales (right panel) . We 
show in addition a parametrization of the total cross section in the left panel. 



factor to be contrasted with approximately energy-independent iT-factor due to next-to-leading 
order corrections. One can see a strong dependence on the factorization and renormalization 
scales which produces almost order-of-magnitude uncertainties and precludes a more precise 
estimation. A better estimate could be done when LHC charm data are published and the 
theoretical distributions are somewhat adjusted to experimental data. 
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Figure 9: Distribution in rapidity (left panel) and transverse momentum (right panel) of c or c quarks at y^s = 

7 TeV. 
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